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ABSTRACT 

Near-infrared interferometers have recently imaged a number of rapidly rotat- 
ing A-type stars, finding levels of gravity darkening inconsistent with theoretical 
expectations. Here, we present new imaging of both a cooler star /3 Cas (F2IV) 
and a hotter one a Leo (B7V) using the CHARA array and the MIRC instru- 
ment at the H band. Adopting a solid-body rotation model with a simple gravity 
darkening prescription, we modeled the stellar geometric properties and surface 
temperature distributions, confirming both stars are rapidly rotating and show 
gravity darkening anomalies. We estimate the masses and ages of these rapid ro- 
tators on L-Rpo\ and HR diagrams constructed for non-rotating stars by tracking 
their non-rotating equivalents. The unexpected fast rotation of the evolved sub- 
giant P Cas offers a unique test of the stellar core-envelope coupling, revealing 
quite efficient coupling over the past ~ 0.5 Gyr. Lastly we summarize all our 
interferometric determinations of the gravity darkening coefficient for rapid rota- 
tors, finding none match the expectations from the widely used von Zeipel gravity 
darkening laws. Since the conditions of the von Zeipel law are known to be vi- 
olated for rapidly rotating stars, we recommend using the empirically-derived /3 
= 0.19 for such stars with radiation-dominated envelopes. Furthermore, we note 
that no paradigm exists for self-consistently modeling heavily gravity-darkened 
stars that show hot radiative poles with cool convective equators. 
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1. Introduction 

Stellar rotation is a fundamental property of stars in addition to the mass and metallicity. 
However it has been generally overlooked in the past century for mainly three reasons. 
Firstly, most stars rotate slowly. Secondly, no complete stellar rotational model is available 
to handle the stellar structure and evolution of a rotating star. Thirdly, although stellar 
rotational velocities in the line of sight vsmi are relatively easy to measure, the inclination 
angles are generally unknown, leaving large uncertainties of stellar rotational velocities. 

While almost all cool stars rotate slowly, rapid rotation is the norm for hot stars. A 
large fraction of hot stars are observed to be rotatin g with equatorial velocities larger than 
120 km s"^ dAbt fc Morrel3 llQQsl : lAbt et al.l[2002l ). Such fast stellar rotation can have 
strong effects on the observed stellar properties. The strong centrifugal forces distort stellar 
shapes and make t hem oblate. Stell ar surface temperatures vary across latitudes due to the 
gravity darkening (jvon Zeipellll924al Jbl). Lower effective gravity at the equator result in lower 
temperatures compared to the poles. This temperature distribution implies that apparent 
luminosities Lapp and apparent effective temperatures T^pp depend on inclination angles, and 
the overall values are hidden from observers. Stellar rotation can also affect the distribution 



of chemical elements, mass loss rate and stellar evolution (IMeynet fc Maederl 12000). Some 
kind of rapidly rotating massive stars may end up as 7-ray bursts (IMacFadyen fc Woosley 
1999h . 



Stellar rotation has been studied mainly through the Doppler broadening line profiles 
in the past, but the obtained information from these studies is limited due to the lack of 
spatial knowledge of stars, such as the inclination angles. An important and reliable way to 
extract such information is through long baseline optical/infrared interferometry, allowing 
us to study the detailed stellar surface properties for the first time. Several rapid rotators 
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These studies have revealed not only the stellar surface geometry but also the surface 
temperature distributions, allowing us to test and constrain stellar models and laws. For 
instance, the surface temperature distributions have confirmed the gravity-darkening law in 
general, but deviate in detail from the standard von Zeipel model (Teff oc where (5 = 
0.25 for fully radiative envelopes). Particularly the studies on Altair and Alder amin prefer 



non-st andard /3 values from the modified von Zeipel model (the /3-free model in IZhao et al. 



( I2OO9I )). These results imply the gravity darkening law is probably only an approximation 
of the surface temperature distribution, the real physics behind is still to be uncovered. 
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In this paper we show our studies of two rapidly rotating stars with extreme spectral 
types in contrast to all the A type stars we have studied: (3 Cassiopeiae and a Leonis, 
observed with the Center for High Angular Resolution Astronomy (CHARA) long baseline 
optical/IR interferometry array and the Michigan Infra- Re d Combiner (MIRC) beain com- 
biner. (3 Cassiopeiae (0 Cas, C aph, HR21) has V = 2.27, jMorel fc MagnenatlllQTsh. H = 
1.584 JCutri et al.lboosh . 1.43 JPucatil l2002l ) . and is located at d = 16 8 pc Jvan Leeuwen 



20071 ). Its mass has been estimated as 2.09 Mrr^ dHolmberg et al.l 120071. see the electronic 



table on VizieR) and it has been classified as F2III-IV ( Rhee et al.l 120071 ). implying it 
was an A type star during main sequence and has evolved - here we will present up- 
dated mass and luminosity estimates (se e Section 5). The rotational velocity has been 

reported between vsin i = 69 km s~^ (iGlebocki fc Stawikowskil |2000| ) and 82 km s~^ 

(IBernacca fc Perinottolll970l ) in the literature, although recent measurernents are more con- 



sistently confined from 69 km s ^ to 71 km s ^ (IGlebocki fc Stawikowskill2000l : lReinersll2006 



Rachford fc Foightll2009l : ISchroder et al.ll2009l ) which we prefer to use in this pape r. Previous 



studies measured its apparent effective temperature range from 6877 Kto 720 K llGray et al. 



200ll : iDaszvhska fc Curieil l2003l : iRhee et al. 



2007; 



Rachford fc Foight 



its radius from 3.43 R^ to 3.69 R^ ljRi chichi fc Percheron 



2002 



2009) and estimated 



Daszyhska fc Cugierl 12003 



Rachford fc Foightl 120091 ). 



a Leoni s (Regulus, H R3982) has V= 1.391 (IKharchenko et a. 



20031), 1.57 flDucatil l2002h . distance d = 24. 31 pc fivan Leeuwen 



rapid ly rotating star, classified as a B7V star (jjohnson &: Morgan 



200% . H= 1.658 fICutri et al 



20071) . It is a well-known 



1953|) or B8 IVn flGrav et al. 



20031 



]). Thet'sin z mea surements from the literature sprea d a large range from ~ 250 km s ^ 



(IStoeckley et al.lll984j ) to ~ 350 km s ^ ( SlettebakI 119631) and w e have adopted here the re- 
cent precise value 317 ± 3 km s~^ from iMcAlister et al.l (120051 ). Regulus is also a famous 



triple star system with the companions B and C forming a b i nary system at ~ 175" away 
from a Leonis A (IMcAlister et al.l 120051 ). Recently iGies et al.l (120081 ) discovered that a Leo- 
nis A is also a spectroscopic binary with a white dwarf company (^ 0.3 M^;^) of the orbita l 
period ~ 40.11 d. The primary mass has been estimated ~ 3.4 M0 (IMcAlister et al.l 120051 ). 
however our study here will show it is much more massive. The diameter of Regulus has been 
estimated several time s in the past because of its brightness and relatively large angular size. 



McAlister et al.l ( 120051 ) combined the CHARA K-band interferometric data and a number of 
constraints from spectroscopy and revealed that Regulus has the polar radius Rpoi = 3.14 ± 
0.06 Roand the equatorial radius Req = 4.16 ± 0.08 R©. 

In this paper, we describe the observations and data reduction in Section 2. Then we 
show the results of /3 Cas and a Leo from both the standard and modified von Zeipel models 
in Section 3. In Section 4, we present aperture synthesis images. In Section 5, the locations 
of the two rapid rotators on L-Rpoi and HR diagrams are discussed. We consider the coupling 
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of the stellar core to the outer envelope and explore gravity darkening from studying these 
two rapid rotators in Section 6. We conclude in Section 7. 



2. Observation And Data Reduction 

The observations were carried out at the Georgia State University (GSU) Center for High 
Angular Resolution Astronomy (CHARA) interferometer array located on Mt. Wilson. The 
CHARA array includes six 1-meter telescopes which are arranged in a Y shape configuration: 
two telescopes in each branch. It can potentially provide 15 baselines simultaneously ranging 
from 34 m to 331 m, possessing the longest baselines in optical/IR of any facility. With 
these baselines, CHARA offers high angular resolution up to ~0.4 mas at the H band and 
~ 0.7 mas at the Khajid. 

The Michigan Infra-Red Combiner (MIRC) is designed to perform true interferometric 
imaging. It is an image plane combiner, combining 4 CHARA telescopes simultaneously to 
provide 6 visibilities, 4 closure phases and 4 triple amplitudes. Currently MIRC is mainly 
used in H band which is further dispersed by a pair of prisms into 8 narrow channels. In 
order to obtain stable measurements of visibility and closure phase, MIRC utilizes single- 
mode fibers to spatially filter out the atmosphere turbulence. The fibers are arranged on 
a v-groove array with a non-redundant pattern so that each fringe has a unique spatial 
frequency signature. The beams exiting the fibers are collimated by a microlens array and 
then focused by a spherical mirror to interfere with each other. Since the interference fringes 
only form in one dimension which is parallel to the v-groove, they are compressed and focused 
by a cylindrical lens in the dimension perpendicular to the v-groove to go through a slit of a 
spectrograph. The data presented here utilized a pair of low spectral re solution prisms with 



R ~ 50. Finally the dispersed fringes are detected by a PICNIC camera (IMonnier et al.ll2004 



20061 ) . The philosophy of the control system and software is to acquire the niaximum data 



readout rates in real time. The details about the software can be found in Pedretti et al. 



(|2009|) 



The integration time is limited by the fast changing turbulence, any turbulence faster 
than 3.5 ms readout speed of the camera will cause decoherence of the fringes. In order to 
calibrate these fringes, calibrators with known sizes adjacent to the targets are observed each 
night. For the acquisition of true visibility, real time flux of the beam from each telescope 



i s also required. Several independent methods (Fiber, Chopper and DAQ , iMonnier et al. 



( 120081 )) are adopted to indirectly measure the 'real time' flux. A recent upgrade of MIRC 
with Photometric Channels has been achieved to directly and more accurately measure the 
flux. Photometric Channels place a beamsplitter right after the microlens array to reflect 
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~ 25% of the flux into multi-mode fibers. The beams exiting the MM fibers go through 
the same doublet and prisms, and hit a different quadrant of the same detector. With 
Photometric Cha nnels MIRC can now measure the visibilities with uncertainty down to 3% 
Jche et allboiol ). 



We observed /3 Cas on 7 nights in 2007 and 2009, and a Leo on 5 nights in 2008. The 
detailed log is presented in Table 1. Figure 1 shows the overall (u,v) basehne coverage of the 
observation of Cas and a Leo . 
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Fig. 1. — Baseline coverage of the all nights observation of /? Cas and a Leo . UV coverage 
can be calculated by dividing these baselines by the wavelength of E band channels. 



Monnier et al. (2007) describes the data reduction pipeline used to process the data, 
which was validated by using data on the calibration binary i Peg. The pipeline first 
computes uncalibrated squared-visibilities and complex triple amplitudes after a series of 
background subtractions, Fourier transformations and foreground subtractions. Then the 
uncalibrated squared-visibilities and complex triple amplitudes are calibrated by the fluxes 
measured simultaneously with fringes. The calibrators with known sizes are observed to 
compensate for the system visibility drift, as listed in Table 2. 



Table 1. Observation logs for /? Gas and a Leo 



Target Obs. Date Telescopes Calibrators 

/3 Cas UT 2007Aug07 S1-E1-W1-W2 7 And 

UT 2007Aug08 S1-E1-W1-W2 a Cyg, 7 And 

UT 2007AuglO S1-E1-W1-W2 ct Cyg, 37 And 

UT 2007Augl3 S1-E1-W1-W2 a Cyg, 7 And, Ups And 

UT 2009Augll S1-E1-W1-W2 7 And, 7 Tri 

UT 2009Augl2 S1-E1-W1-W2 7 And, 7 Tri 

UT 2009Oct22 S2-E1-W1-W2 37 And, v And, e Cas, r? Aur 



a Leo UT 2008Dec03 S1-E1-W1-W2 d Leo 

UT 2008Dec04 S1-E1-W1-W2 54 Gem, t] Leo 

UT 2008Dec05 S1-E1-W1-W2 B Hya, 9 Leo 

UT 2008Dec06 S1-E1-W1-W2 54 Gem, B Hya, »7 Leo 

UT 2008Dec08 S1-E1-W1-W2 B Leo 



Table 2. Calibrator diameters 



Calibrator 


UD diameter ( mas) 


Reference 


7 And 


0.659 ± 0.017 


b, c, d 


37 And 


0.682 ± 0.030 


b, c 


u And 


1.14 ± 0.007 


a, b, c, d 


o- Cyg 


0.542 ± 0.021 


a 


7 Tri 


0.520 ± 0.0125 


b 


e Cas 


0.351 ± 0.024 


c, d 


r? Aur 


0.419 ± 0.063 


c 


6» Leo 


0.678 ± 0.062 


b, c 


r] Leo 


0.644 ± 0.068 


c 


54 Gem 


0.735 ± 0.033 


b, c 


6 Hya 


0.463 ± 0.031 


c, d 



HMerandl ||2008|') 
feervella fc Fouguel ||2008|'| 
iBarnes et ai] j 19781) 
' ^Bonneau et al.l 1I2OO6I') 
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3. Modeling of Rapid Rotators 

We construct a 2D stellar surface model in this paper: the modified von Zeipel model. 
The model contains six free parameters, stellar polar radius, the polar temperature, the 
ratio of angular velocity to critical speed oj / Wcrit, the gravity darkening coefficient the 
inclination angle, and the position angle (east of north) of the pole, to describe the stellar 
radius, surface effective gravity and temperature distributions across stellar surface. The 
mass of a star is given and ffxed in each model fftting process. Given the stellar mass, stellar 
polar radius an d oj / ujrr\t, the stellar ra dius and surface effective gravity at each latitude can 



be determined (lAufdenberg et al.ll2006l ). Then given the stellar polar temperature and /3, the 
stellar surface temperature distribution can be computed from the gravity darkening law(T oc 
^fj). Lastly the orientation of the star is described by the inclination angle and position angle. 
In the model, we assume the solid-body rotation for simplicity; a more complicated and 
realistic model would consider the differential rotation which requires additional information 
(such as spectral lines) for fitting. The gravity darkening coefficient /3 is a free parameter 
in the model. By fixing /?, the model reduces to the standard von Zeipel model {j3 = 0.25, 
radiative case) or Lucy model (/? = 0.08, convective case). 



In earlier work (IMonnier et al.l 120071 ) , we found that allow /3 to be a free parameter 
greatly improved the fit to the interferometric data. This fiexibility allows us to indepen- 
dently test the validity of the standard von Zeipel and Lucy prescriptions. Furthermore, the 
mixture of radiative and convective regions in the same star may also cause deviations from 
expected values. For example, the polar temperature could be thousands of degrees higher 
than the equator temperature, resulting in a situation where upper atmosphere may be ra- 
diative at the poles while convective at the equator. In general, the value of also depend s 



on various approximations made for the atmosphere, radiation transfer etc. ( lClaretlll998l ). 
Therefore in our modified von Zeipel model, instead of setting (3 to be fixed, we allow (3 to 
change as a single free parameter of the model to fit the interferometric data. For compari- 
son, we also present models with (3 fixed to the appropriate standard value. The error bars 
of stellar parameters from the modified von Zeipel model are in general larger than those 
from standard von Zeipel model or Lucy model. This is because there are certain degrees 
of degeneracies between the gravity darkening coefficient /3 and other stellar parameters, as 
discussed below. 

During the model fitting process, the modified von Zeipel model is converted into a 
projected stellar surface brightness model, which is constrained by the observed V and i/0 



^We used H magnitudes and errors from only 2MASS catalogue to constrain the model fit ting. After 



we submitted the paper, we found more precise measurements of H magnitudes (jDucatil 120021 ) which are 
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band photometric fluxes and three kinds of interferometric data from each night: squared 
visibihties, closure phases and triple amplitudes. In the modified von Zeipel model, the stellar 
surface is divided int o small patches. The intensity of each patch is computed from Kurucz 
model (Kurucz 1992 1^ given the temperature, gravity, viewing angle and wavelength, so that 



the modified von Zeipel model can be converted into the projected brightness model. The 
projected brightness model is then converted into the same three kinds of interferometric 
data above by direct Fourier transform to fit to the observed data. We use 4 sub-bands 
(binning two adjacent narrow channels dispersed by the MIRC prisms) across H band for 
accuracy. In addition, the apparent Fand if band photometric fluxes are obtained from the 
projected brightness model to fit to the observed values. Observed v sin i is not directly used 
in the model fitting, but used to cross-c heck the results from model fitting. The detailed 



process is described in IZhao et al.l ( 120091 ) and reference therein 



Data errors consist of random errors, errors due to variation of seeing condition, and 
calibration errors from using incorrect diameters of the calibrator targets. To get the errors 
from the first two parts, we treat the data from each night as a whole package and bootstrap 
packages randomly with replacement. Then we fit the sampled data and repeat fifty times 
to get the distribution of each model parameter. The upper and lower error bars quoted 
here are such that the interval contains 68.3% probability and the probability above and 
below the interval are equal. For the error from the third part, we used simple Monte Carlo 
sampling using the our estimated angular size uncertainties - these errors turned out to be 
somewhat smaller than the error from the first two parts. 

We should point out that the stellar mass has to be given and fixed at the beginning 
of each model fitting process, but at first does not agree in detail with the model estimated 
from the fitting results on both L-R^^x and HR diagrams using the rotational correction (see 
Section 5). Our approach here has been to adopt the mass from the literature for the first 
attempt in the model fitting. The mass estimation from the first attempt is then used in 
the second round of model fitting process. This procedure is repeated until the mass given 
in the model agrees with what comes out of the model fitting. The final mass is referred as 
the model mass in our paper. The stellar metallicities are adopted from the literature and 
fixed throughout. The distances of the targets are also adopted from the literature. 

We also calculate the stellar mass based on the measured v sin i ra nge from the litera- 



ture, which is referred as the oblateness mass and was first proposed by IZhao et al.l ( l2009l ). 



For each bootstrap, we extract the inclination angle, polar radius and uj / Wcrit from the 



consistent with our model values within l-cr. 
^Data downloaded from kurucz.harvard.edu/ 
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best fitting, then uniformly sample vsini values 100 times in tlie given range to obtain a 
mass distribution. By combining tfie mass distribution from each bootstrap, we obtain the 
whole mass distribution from which the upper and lower mass bound can be calculated such 
that the interval contains 68.3% probability and the probability above the upper bound and 
below the lower bound are the same. To compute the best estimation of the stellar mass, we 
use the best estimations of the inclination angle, polar radius and u / cjcrit from the model 
fitting of all nights, and the v sin i value to be the mean of the measured range from the 
literature. 



3.1. /3 Cassiopeiae 



We ad opted the following basic properties of /3 Cas from the literature as input s: distance 
= 16.8 pc Jvan LeeuwenI boOTh and metallicity [Fe/H] = 0.03 flCrav et al.lboOlh . We take 
[Fe/H] = which is t he closest value to the observation to extract intensities from Kurucz 
model. M = 2.09 (IHolmberg et al.l 120071 . see the electronic table on VizieR) is adopted 



for the first attempt of the model fitting. The fitting results and final parameters from the 
modified von Zeipel model are shown in Fig. 10 in the Appendix and the middle column 
of Table 3 respectively. The results show that (3 Cas is rotating more than 90% of its 
critical rate, which causes its radius ~ 24% longer at the equator than at the poles. The 
temperature at the pole is about 1000 K higher than that at the equator. These significant 
differences between the poles and equator imply that the Lapp and 7|pp are highly dependent 
on viewing angles. The best model mass estimation of its non-ro tating equivalent from L-R^ca 
and HR diagrams is 1.91M0(Fig. 6), lower than 2.O9M0from lHolmberg et al.l (120071 ). The 
oblateness mass estimation from vsmi range 69 km s^^ to 71 km s~^ is 1.771*^005 M©, 
which is consistent with our model mass within the error bars. j3 = 0.146 from the modified 
von Zeipel model fitting is significantly different from standard values for either radiation- 
dominated or convection-dominated envelopes. The inclination angle is low, implying we are 
looking at more the polar area than the equatorial shown in Figure 4 (see Section 

4). This is why the apparent luminosity Lapp is higher than Lboi • 



Claretl (120001 ) has computed the evolution of gravity darkening coefficients for different 



stellar masses, and showed that at such low T^s as f3 Cas it should be convection-dominated 
in the envelope. Fixing gravity darkening coefficient /? = 0.08 (Lucy model) for convective 
envelopes, we run model fitting again and the results are shown in the right column of 
Table 3. The best fitting x^s for this model is much worse, nearly a factor of 2 higher. 
Many parameters from the Lucy model are similar to those from the modified von Zeipel 
model, except the temperature at the equator. This is not surprising because the low /? 
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value means the weak dependence of the temperature on gravity, namely the temperature 
at the equator will be closer to that at the poles for the Lucy model. Consequently the 
luminosities and temperature Lapp , Z/boi and T^pp are a httle higher than those from the 
modified von Zeipel model. The modified von Zeipel model gives significantly lower than 
the Lucy model, especially that from the closure phase data which is sensitive to asymmetric 
structures on the stellar surface. This implies the modified von Zeipel model describes the 
surface temperature distribution better, ruhng out the Lucy model in this case. This is 
also confirmed by comparing the model vsini with the observed values: vsini — 72.413 5 
km s~^ from the modified von Zeipel model agrees with the observation 69 km s~^ to 71 
km , while from the Lucy model vsini = 81.3l?;o km s~^ deviates strongly from the 
observation. Further more, the oblateness mass and model mass don't agree with each other, 
suggesting that the Lucy model is not self-consistent in this case. 

We found that the low inclination angle induces strong degeneracies between some 
parameters during the model fitting. For example when a star is pole-on the darkness 
at the equator could be due to either the high angular velocity or the high gravitational 
darkening coefficient since the oblateness can not be directly constrained from this viewing 
angle. Therefore we explore the probability spaces of gravity darkening coefficients j3 with 
inclination angles and oj / ojcnt to assess possible correlations. For example, we first search 
the best model fitting results of all nights on a 40 x 40 grid of (5 and inclinations by fixing 
these two parameters on each pixel. Generally if these two parameters are independent, 
then the probability of their true values falling into each pixel is oc e~°'^^^ . However in this 
case the two parameters are dependent, we modify the probability oc e~°^ , where o; is a 
variable to be determined. Then we overplot the results of the two parameters from each 
bootstrap onto the probability space (not shown in the figure), and find the contour of the 
same containing 68.3% of bootstrap results, from which a can be computed. The contour 
is defined as l-cr. 

The left panel of Fig. 2 shows the degeneracy between /3 and the inclination. The 
contour represents the 68.3% probability level, and is weakly elongated in one direction. We 
further overplot onto the probability space the observed v sin i range which intersects the 
contour. This means a precise v sin i measurement would significantly constrain the stellar 
parameters from our model fitting. The same idea is appUed to the probability space of p and 
I ^^crit (Figure 2 right) which shows a stronger correlation between these two parameters. 
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Table 3. Best-fit and physical parameters of /5 Cas 



Model Parameters Modified von Zeipel model (/3-free) Lucy model (/3 = 0.08) 

Inclination (degs) IS-St}'! '^^■^-o'.l 

Position Angle (degs) -7-09^2 40 -l-Sl? ? 

Tpoi (K) 7208lg TlOStit 

Rpol ( mas) 0.849t° °^3 0.8351°;°?^ 

Lo I u;,Ht 0.920t0;°24 CQSO^jOIqIso 

0.146t°;o" 0.08 (fixed) 

Derived Physical Parameters 



Tcq (K) 


6167+2? 


6487l}7 


Rcq ( R©) 


70+0.10 

'^-0.09 


o 77+0.16 
'-0.04 


Rpol ( R©) 




r,i+0.13 
'-'•'^J^-0.04 


Bolometric luminosity Lbol (Lg) 


^-'-■■'-0.7 


' -0.3 


Apparent effective temperature T^pp ( K) 


6825 


6897 


Apparent luminosity Lapp (L©) 


27.3 


28.3 


Model vsvai ( km s"-"^ )^ 


72.4+3^ 


79.8+?-^ 


Rotation rate (rot/day) 


1 1 9+0.03 
^■^^-0.04 


1 16+°"1 
J^-J^O-0.06 


Model mass (M©)'' 


1.91 ±0.02 


1.95 ±0.03 


Oblateness mass (M©) 


T 77+0.17 
^- ' '-0.05 


1 4Ci+0.12 
-'-■^'^-0.27 


Age (Gyrs)*" 


1.18 ±0.05 


1.09 ±0.03 


Model y Magnitude'^ 


n 904+0.012 
^•^°*-0.019 


9 9cr-, +0.020 


Model H Magnitude'^ 


, .,00+0.007 

1.0»0_Q ggy 


1 304+0.010 


of various data 


Total xi 


1.36 


2.53 


Vis2 xl 


1.26 


1.56 


CP xl 


2.18 


4.81 


T3amp xl 


0.45 


0.60 


Physical Parameters from the literature 


[Fe/H]<= 


0.03 




Distance (pc)' 


16.8 





''Observed D sin i =69 km s ^ to 71 km s 1 ijGIebocki &: Stawikowskill200Cl:lReinerj20od : lRachford fc Foighd 
I2OO9I : ISchroder et aLllSooj) 

''Based on the stellar evolution model l iYi et al.ll200ll . bOOSi : iDemargue et alj|2004h 

iZhao et al.) l|2009h 

I Morel fc Magnenaj|l978l , with arbitrary error), Hmag = 1.584 ± 0.174 iCutri et all 



y mag = 2.27 ± 0.01, I Morel 
l2003f) , 1.43 ± 0.05 llDucatill20o3i ^ 



^Grav et al.] l l200ll 'l 
Ivan LeeuwenI ||2007| ') 
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3.2. a Leonis 



We first fit the stellar surface of the modified von Zeipel model to the interferometric 
data of a Leo . Th e parameters we a dopted from the literature are given as fol lowing: 
distance = 24.31 pc dvan Leeuwenl boOTh. metalhcity [Fe/H] = 0.0 JCrav et al.lboosh . Mass 
= 3.4 ( iMcAlister et al.ll2005l ) was used for the first attempt of the model fitting. The 
fitting results from the modified von Zeipel model are shown in Fig. 11 in the appendix, with 
the final stellar parameters listed in the middle column of Table 4. a Leo is rotating at 96% 
of its critical speed, causing the equatorial radius about 30% longer than the polar radius. 
The temperatures at the poles are more than 3000K hotter than that at the equator. The 
gravitational darkening coefficient /5 from the fitting is again different from the "standard" 
values for either radiative or convective envelopes. The results show that a Leo is almost 
equator-on, which is shown as a dark strip in Figure 5 (see Section 4). Therefore the Lboi is 
higher than the Lapp . The model mass from HR diagram is 4. 1 5 ± .06 Mq. Adopting the 
vsini range vsini = 317 ± 3 km s^^ from iMcAlister et al.l (|2005[ ) paper, the oblateness 



mass estimation corresponding to the model mass is 3.661q'28M0, which also agrees with 
the model mass within the errors. The large errors of the oblateness ma ss is due to the 



degen eracy of stellar parameters as discussed later. The observed vsini (IMcAlister et al. 
20051 ) is consistent with our derived value 336^24 km s~^ with error bars. 



Theoretically the high surface temperature of a Leo suggests that the envelope is fully 
radiative, corresponding to the gravity darkening coefficient /3 = 0.25. We fit the model 
again using the fixed /3 value, which is the standard von Zeipel model. The best fitting x^s 
for this model is much worse, nearly a factor of 2 higher. For completeness, we have included 
the results in the right column of Table 4. In this scenario, a Leo is rotating even faster. The 
larger gravitational darkening coefficient and faster rotation imply even larger temperature 
difference between the poles and equator. However the derived equatorial temperatures from 
the modified and standard von Zeipel models agree with each other. This is because Regulus 
is almost equator-on, the observed interferometric data is dominated by information from 
the equator. The x^s of the various interferometric data from the modified von Zeipel model 
are all significantly smaller than those from the standard von Zeipel model, supporting the 
modified von Zeipel model with /? = 0.19 is preferred to describe the surface properties of 
Regulus, ruling out the standard von Zeipel value. This conclusion is also supported by the 
disagreements between the model mass and oblateness mass from the standard von Zeipel 
model, and between the model and observed vsini values (see Table 4). 

We expect some degeneracies of parameters from the modified von Zeipel model fitting 
because of the symmetry of the equator-on orientation. Two figures of probability space of a; / 
ojcritand the inclination vs. f3 are shown in Figure 3. Both pictures show a strongly elongated 
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contour of the probability, implying significant correlation between these parameters. The 
solid contours s h ow th e 68.3% probability. We overplot the observed vsini range from 



McAlister et al.l (120051 ) . which intersects the contour with a much smaller common area. 



Therefore a precise v sin i measurement would significantly reduce the degeneracy between 
the parameters and constrain them much better. 



Based on only visibility data, iMcAlister et al.l ( l2005l ) modeled a Leo and our new model 
results are generally consistent with this earlier work. Since MIRC has higher angular res- 
olution, better UV coverage and the closure phase data, our data is more sensitive to the 
detailed structures such as the inclination and position angles. This work found acceptable 
fits for (3 values between 0.12 and 0.34 (best fit at 0.25), a range consistent with our more 
refined analysis. Although our estimations of the bolometric luminosity Lboi of Regulus is 
similar to those from their paper, the HR diagram (Fig. 7) from our results suggests that 
the mass of the non-rotatir ig equivalent of Regulu s is 4.15 ± 0.06 M©, much more massive 
then the 3.4 ± O.2M0that iMcAlister et al.l (120051 ) obtained using the surface gravity log g 
from spectral analysis. Their results show that the non-rotating equivalent of Regulus has 
lower m ass and conseque ntly lower Lboi than rapidly rotating Regulus, which is in contrast 
to what ISackmannI (Il970l ) found that a non-rotating equivalent actually has higher Lboi than 
its rapidly rotating equivalent. 



4. Imaging Of Rapid Rotators 

Interferometric data contains information of the Fourier Transform of the projected 
surface brightness of sources. Therefore, in theory, a stellar image can be reconstructed 
directly from the data. But in reality because of the finiteness of UV coverage and uncertainty 
from measurements, many different images fit well to the same interferometric data. We use 



the ap plication 'Markov-Chain Imager for Optical Interferometry' (MACIM; llreland et al. 
(120061 )) to construct images for /3 Gas and a Leo . It is usually difficult to image nearly 
point-symmetric objects because the closure phases will be close to either or 180 degrees, 
making it harder to constrain the detailed structure. /? Gas is close to pole-on and a Leo is 
almost equator-on, which are two cases of the point-symmetry. 

One strategy to image these kinds of stars is to take advantage of some prior knowledge. 
Stars are confined in certain area with elliptical shapes approximately. Therefore we employ 
a prior image which is an ellipse with uniform surface brightness. The spatial and geometric 
parameters of t he ell ipse come from the model fitting. The detailed process can be found in 



Monnier et al.l (120071 ). 
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Table 4. Best-fit and physical parameters of a Leo 



Model Parameters 



Modified von Zeipel model (/3-free) von Zeipel model (/3 = 0.25) 



Inclination (dogs) 
Position Angle (dcgs) 
Tpoi (K) 
Rpoi (mas) 

/3 



Derived Physical Parameters 



86.3 
258 
14520 
0.617 
0.962 
0.188 



+1.0 
-1.6 
+2 
1 

+550 
690 
+0.010 
0.009 
+0.014 
0.026 
+0.012 
0.029 



1+0.07 
-0.06 



87.5"* 



c+0.2 
-0.1 

2591^ 
16190liio 
0.605+?,°oj 



0.001 
-0.001 
-0.002 



0.969 
0.25 (fixed) 



10920+ 

4 ,-+0.007 
^•^'-0.006 

q if; + 0.005 

•^•J^ '-'-0.004 
431^1* 
12650 
294 
346t^ 

1 -p+O.Ol 

^- "J-o.oi 
4.52 ±0.05 
o 44+0.08 
''■^^-O.Ol 

0.05 ±0.01 

-, 020+0.017 
i.ozy_Q Q21 

l.OOU_0^015 



Tcq (K) 
Roq ( R© ) 

Rpol (R©) 

Bolometric luminosity Lbol (L©) 

Apparent effective temperature T^pp ( K) 

Apparent luminosity Lapp (Lg) 

Model tjsini ( km 

Rotation rate (rot/day) 

Model mass (M©)'' 

Oblateness mass (M©) 

Age (Gyr)b 

Model y Magnitude'* 

Model E Magnitude'* 



4.21"' 



+0.05 
-0.04 
+27 
28 



3.22 
34i:l 
12080 

252 
3361^^ 

J^-°*-0.04 

4.15 ±0.06 

3 fifi+O-79 

0.09 ±0.02 

1 000+0.002 

1.0yO_Q ggg 

1 c:7o+0.004 



of various data 


Total xi 


1.32 




2.57 


Vis2 xl 


0.76 




1.26 


CP xl 


1.97 




3.80 


T3amp xl 


0.92 




1.52 


Physical Parameters from the literature 


[Fe/H]'= 




0.0 




Distance (pc)* 




24.31 





''Observed vsmi = 317 ± 3 km s"! l|McAlister et al.ll2005h 

''Based on the stellar evolution model liYi et al.ll20"0ll . |2003|: iDemargue et al.ll2004h . 
' IZhao et al.) bOO^ 

** V "mag = 1.391 ± 0.007 llKharchenko et aLllioOol) . i/mag = 1.658 ± 0.186 jCutri et al.ll2003l) . 1.57 ± 0.02 llDucatil 
l2002h 

■ iGrav et al.] l l200ll'l 

Ivan Leeuwenl ||2007| ') 
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The left panel of Figure 4 shows the reconstructed image of {3 Cas . The reduced of 
the image is 1.20, comparable to our best-fit models. We overplot longitudes and latitudes 
with solid lines from the model and include contours of surface brightness temperatures with 
dashed lines. The right panel shows the image from the model fitting, overplotted with the 
surface brightness temperature contours from the model. Because of the inclination angle, 
the surface brightness temperature contours do not coincide with latitude contours. We find 
that the two images are consistent with each other in general. The images show a center 
bright region which is one pole of /3 Cas . The surface brightness drops gradually towards 
the edge due to gravity darkening. One may also notice limb-darkening at the edge of the 
stellar image. 

The left panel of Figure 5 shows the image of a Leo with latitudes and longitudes from 
the model, and surface brightness temperature contours. The reduced of the image is 0.78. 
The right one shows the image from model fitting. As opposed to (3 Cas , a Leo is almost 
equator-on and the dark equator stretches along the North-South direction. One noticeable 
phenomenon is that the poles are not located exactly in the hot region. This is because in 
this particular case the poles at the stellar image edge look cooler due to lime-darkening, 
causing the brightest regions to shift towards the center of the image. 



5. Stellar Evolution Tracks 

One interesting topic for rapidly rotating stars is to locate their positions on the Hertzsprung 
-Russell (HR) diagram and compare with stellar models. This topic contains two issues. 
First, traditional photometry observations only see the apparent luminosities Lapp and ap- 
parent effective temperatures which depend on stellar inclination angles; the bolometric 
luminosities Lboi of rapid rotators are hidden from the observers. Interferomet ric observa- 



tions allow us to construct 2-D surface models of stars, thus to obtain the Lboi (IZhao et al. 



20091 ) . We obtain the gravity and temperature distributions across the stellar surface from 
the model fitting. From Kurucz model, we are able to retrieve intensities from each patch 
of stellar surface, and then integrate the radiation all over the star to obtain the bolometric 
luminosity Lboi □• By comparison we also compute an inclination curve which shows stellar 
Lapp and T^pp as a function of the inclination angle, and we can mark the one corresponding 
to its inclination from the model fitting. The Lapp can be calculated by Lapp = 4:HcPF]^oi, 



•^The "overall effective temperature" T^fi can be estimated from the iboi divided by the total surface 
area; However, in the case of a rapid rotator, this overall effective temperature is just a definition with 
limited physical meaning, so it is not used to infer the masses or ages of stars in this paper. 
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where d is the distance and Fboi is the bolometric flux computed by integrating flux from 
each grid over the projected area. Then the T^pp is obtained by cr(T^pp )^ = yrrf^Fboi/^p 
where ^proj is the projected area. 



^proj ; 



Second, typical HR diagrams are constructed for non-rotating stars, it is inappropriate 
to place a rapid rotating star on such diagrams. A rapidly rotating star shows a little lower 
Lboi than L^j. from its non-rotating equivalent (an imaginary spherical star which a rapid 
rotator would turn out to be if it spins down to no angular velocity), meaning a rotating 
star will evolve as a lower mass star on HR diagram. Therefore the interpreted mass and 
age from the rotating star deviates from the true values. To partially solve this problem, 
one has to convert the properties of a rapidly rotating star to its non- rotating equivalent. 
Studies have shown that the bolometric luminosity and polar radius do not change much 
as a star spins up. Following this, we alter the traditional HR diagram to a new one with 
axes of bolometric luminosity and polar radiu s (L-Rpr,] diagram), a nd locate rotating stars 
on the new diagram to infer the mass and age (jPeterson et al.ll2006l . private communication, 
2010). To compare with the astronomy-friendly HR diagrams, one can also translate these 
two values of non-rotating equivalents into L^^ and . 

The left panels o f Fig. 6 and 7 show C as and a Leo on L-R^ox diagrams from Y'^ model 



(lYi et al.ll200ll . l2003l : iDemarque et al.ll2004j ). The cross and square symbols represent the 



bolometric lumiii osity and polar radius before and after the rotational correction respectively 
(ISackmannI Il970l ) . The corrections are trivial: L^^ and Rpoi,nr decrease by 5.5% and 1.3% 
respectively for a 2 solar mass star as it spins up to close to critical speed. So on L-R^oi 
diagrams one may even directly use Lboi and Rpoi of a rotating star for rough interpretations 
of its mass and age. We have begun work on a more exact formulation using a new grid of 
rotating models, but this is the subject of future detailed paper. 

The traditional HR diagrams are shown in the right panels of Fig. 6 and 7. The solid 
lines are the inclination curves, which show the Lapp and T^pp as a function of inclination 
angles. The star symbols on the curve represent the estimated inclination angles. The 
square symbols stand for L^^ and 1^ of the non-rotating equivalent. The position of non- 
rotating equivalent on HR diagram deviates severely from the position of the rapidly rotating 
equivalent based on its apparent values. For instance, Regulus would be about 0.08 Gyr 
older and 0.5 MqIcss massive from its Lapp and T^pp than from Lnr and . So we strongly 
recomme nd to correct for th e effects of rotation when placing a rapidly rotating star on HR 
diagram. IZhao et al.l (120091 ) didn't adopt this correction, which may lead to an additional 
error in determining age and mass of rapidly rotating stars. 
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6. Discussion 

6.1. Stellar Core-Envelope Coupling 

Measuring uj / oJcrit as a function of age provides a way of studying the coupling between 
the stellar core and envelope in terms of angular momentum. As a star evolves along the 
main sequence, the core contracts and spins up due to the conservation of the angular 
momentum, while the spherical-shell envelope expands and spins down, ujcnt also drops as 
the star expands. Given the initial rotational conditions and the evolution of stellar inner 
structure, the evolution of a; / Wcrit depends only on how much the core and envelope are 
coupled. In the case when the core and envelope are not coupled, the angular velocity of the 
envelope changes roughly proportional to R~^. The ratio oJcrit is proportional to R~^-^. So 
UJ I oJcrit decreases roughly as R~°'^ as a star expands. While in the other extreme case of 
solid body rotation, namely the core and envelope are fully coupled, the core transfers the 
most angular momentum to the envelope, and uj / Ucrix, may increase as a star expands. We 
can also predict its value in the past, knowing the current u / Wcrit- 

One critical component in the discussion above is the evolution model of stellar inner 
structure. While several such models are available for non-rotating stars, we can not find one 
for the general case of rotating stars. We justify that a non-rotating stellar model is a good 
approximation for calculating evolution of internal density profiles because rotation has very 
little effect on iso-potential surfaces inside the star. For instance, a rapidly rotating star with 
^ I i^crit = 0.9, its equatorial radius is elongated by only 21.6%, but gravity quickly dominates 
as one looks deep into the star. This means Wcrit is much larger than angular velocity at 
certain radius and smaller, and the structure can again be approximately described by a 
non-rotating stellar model. So in the following calculation we adopt a non-rotating stellar 
mode0 

By computing how the moment of inertia changes with time, we are able to calculate the 
evolution of a; / Wcritfor a 1.9 M0 non-rotating star (Fig. 8). In the left panel, all the values 
are normalized to their initial values. The solid line shows the evolution of the stellar radius, 
the dotted and dashed lines show the evolution of the ratio uj / Ucrit when the core and 
envelope are fully coupled and uncoupled. When the core and envelope are uncoupled, the 
ratio drops as the star expands as expected. When the core and envelope are fully coupled, 
the ratio actually increases a little due to the transference of angular momentum from the 
core to the envelope. This result may explain high uj / Wcrit value of /3 Cas . 



"'EZ-Web http : / / www, astro . wise . edu / ^townsend/st atic . php?ref =ez-web is a web-browser interface to the 
EZ evolution code (IPaxton 2004 ). developed and maintained by Rich Townsend. 
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In the right panel, we use the ratio uj / Wcrit = 0.92 from model fitting as the current 
value of [3 Cas , and trace back to its previous values in the extreme cases of full-coupling 
and no coupling. We notice that if the core and envelope are not well-coupled (dashed line), 
the ratio will exceed the unit in the past, which is not allowed. On the other hand if they 
are totally coupled (dotted line), the ratio value remains below 1. Reading off the panel, u 
I Wcrit changes more rapidly in the past ~ 0.5 Gyr if the core and envelope are not coupled. 
These results suggest that during the stellar evolution of /3 Cas , the angular momentum is 
efficiently transferred from t he core to the envelope iri the past 500 Myr. These results seem 
to confirm earlier findings by lDanziger &: Faberl (119721 ) based on analysis of v sin i statistics. 



6.2. Gravity Darkening CoefRcient 

Von Zeipel brought up the idea of gravity darkening in 1924 and predicted the standard 
value of (3 to be 0.25 for stars with fully radiative envelope. Our group have studied five rapid 
rotators (a Aql, a Cep, a Oph, a Leo, (3 Cas) up to now, four of them show non-standard 
Gravity darkening coefficient (/?) values from the modified von Zeipel model fitting, a Oph 
was only fitted with /3-fixed model because of the high degeneracy between gravity darkeniri g 



coefficient and rotational speed due to its almost equator-on orientation (jZhao et al.ll2009l ). 



In Fig. 9 we plot the results of /3 versus temperature for the four targets with their grav- 
ity darkening coefficients obtained from the modified von Zeipel model fitting. The shadow 
areas show the temperature ranges from the pole to equator and the l-cr uncertainties of /5 
from the model fitting for each star. For comparison, we also plo t the so l id line representing 



the predicted relation between /3 and temp erature adopt ed from IClaretl (120001 ). We digitize 



the evolution plot of a 2 solar mass star in I Clare tl (120001 ) paper and extend /3 to high tem 



perature 14500 K with /5 fixed to 0.25. We should point out that the predicted relation shifts 
a little to lower temperature for stars with higher masses, but it is not a big issue in our 
case. For a Cep, a Aql and (3 Cas, their masses are close to 2 Mq, so they can share the 
same relation, a Leo is much more massive than 2 Mq, the predicted curve shifts to low 
temperature a little (less than lOOOK). 

Fig. 9 shows that a Cep, a Aql and /? Cas partially intercept the transition area of the 
predicted curve, meaning that the equatorial regions might star to show convection. In our 
model fitting, we use a single f3 to describe the relation between the gravity and temperature, 
instead of letting f3 change as a function of temperature. This may partially explain why these 
three stars have non-standard /? values, because their poles could be radiation-dominated 
while the equators convection-dominated, the resulting (3 may be some weighted values across 
the stellar surfaces. However the analysis here is non-physical, a detailed stellar model that 
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includes radiation and convection in a rapidly rotating star is required to fully understand 
the gravity darkening law of these stars with intermediate temperatures. 

However a Leo has such high temperature range that even the equator is supposed to be 
fully radiative theoretically. So the poles and equator will share the same = 0.25, justifying 
the standard von Zeipel model in this case. But our result still prefers non-standard /3 = 
0.1881q;q29 • O'^^ pos sible exp l anatio n is that even at such high temperature, the envelope 



is not fully radiative. iTassoull ( 120001 ) concludes that sohd-body rotation is impossible for a 
pseudo-barotrope in static radiative equilibrium. The sohd-body rotation will disrupt the 
constancy of the temperature and pressure over the stellar surface, and cause the temperature 
and pressure gradients between the equator and poles. The gradients will induce a flow of 
matter which forms a permanent meridional circulation and break down the strict radiative 
equilibrium. The matter flow may further lead to the failure of our model assumption: 
solid-body rotation. The material from higher latitudes carries less angular momentum than 
those from lower latitudes. The meridional flows moving towards higher or lower latitudes 
will speed up or slow down the rotational speed of local material on their way, which triggers 
differential rotation. 



Another study from iLovekin et al.l ( 120061 ) compares the effective temperature distribu- 



tion across the surface of a 6.5 Mq solid-body rotator between a stellar evolution model with 
rotation (ROTORC) and von Zeipel's law, and finds that the temperature distribution is 
shallower in the model which is co nsistent with lower /j value we o btained from a Leo . A few 



observations on W UMa systems (IKitamura &: Nakamuralll988l : iPantazis fc Niarchoslll998[ ) 



roughly confirm von Zeipel's law, but with very large scatter. The material flows on the sur- 
faces of these stars are less complicated due to an important feature of the binary systems: 
the stars are tidally locked by their companions. Hence the stellar differential rotations are 
effectively depressed and the resulting solid-body rotations are well regulated. Therefore 
these stars may maintain radiation-dominated envelopes which validate the standard von 
Zeipel model. 

Based on the similar /3 value found for all our objects and for a Leo in particular, we 
recommend researchers adopt a new standard /3=0.19 for future modeling of rapid rotating 
stars with radiative envelopes. 



7. Conclusion 



We have studied two rapid rotators with extreme spectral type: /3 Cas and a Leo 
observed by CHARA-MIRC. By fitting the modified von Zeipel model, namely the solid- 
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body rotation model with free-/? gravity darkening law, to observed infrared interferometry 
data and V and H photometric fluxes, we find both stars are rotating at close to critical 
speed: u / Wcrit= 0.92 and 0.96. The fast rotations elongate their equators by 24% and 30% 
compared with their poles, and their equatorial temperatures are lOOOK and 3000K cooler 
than their polar values. We estimated the mass of a Leo to be 4.15 ± 0.06 M0 from both 
L-Rpoi and H R diagrams corrected f or rotational effect, and it is much higher than 3.4 ± 0.2 



M0found by lMcAlister et al.l (120051 ). We have also reconstructed aperture synthesis images 
using MACIM. The images are consistent with the temperature distribution from the model 
fitting. 

We discussed the evolution of u / Ucrit- The ratio could increase or decrease depending 
on how much stellar cores and envelopes are coupled. In the case of fully coupling, u / oJcrit 
increases a little during main sequence and sub-giant branch due to the angular momentum 
transferred from the core to the envelope. Our study on /3 Gas, which is about 1.18 Gyrold 
but still rotating at 92% of its critical speed, suggests the core and envelope are well coupled 
during the evolution. 

All our targets from the modified von Zeipel model fitting prefer the non-standard grav- 
ity darkening coefficients, especially in the case of a Leo whose envelope should be fully 
radiative because of the high surface temperature range IIOIOK - 14520K. One possible 
reason is that solid-body rotation breaks down the constancy of temperature and pressure 
on the stellar surface and induces meridional flow, which violates strict radiative equilib- 
rium. Furthermore the meridional flow may result in differential rotation which causes the 
failure of our solid-body rotation assumption. To explore this possibility in the future, we 
will construct a differential rotation model to fit observed high resolution spectra of these 
rapidly rotating stars. Until better models are created, we recommend using the empirically- 
determined gravity- darkening coefficient (3 = 0.19 for rapidly-rotating stars with radiative 
envelopes. 
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A. Appendix 

Visibility data and fitting results from one single night arc shown in Fig. 10 and 11. 
The upper left panels show all 7 and 5 nights combined visibility data of j3 Cas and a Leo re- 
spectively, overplotted with the visibility curves of uniform disks with diameters of major 
and minor axiscs from model fitting. The other panels show the model fitting and imaging 
results compared with one single night data. The date of that night of /3 Cas is Oct. 22nd 
2009, and that of a Leo is Dec 8th 2008. 
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Fig. 2. — Probability spaces of /? Cas show the degeneracy between stellar parameters. 
The left panel shows the probability space of the gravity darkening coefficient /? and the 
inclination angle; the right one shows that of /3 and the fraction of critical angular velocity 
u I ojcrit- The solid contours represent the l-a levels, containing 68.3% of the probability. 
And the dashed lines connect pixels in the probability space with the same vs\.ni values 
from model fitting. The value range 69 km s~^ to 71 km s~^ is adopted from the literature, 
and the corresponding lines intersect the 1-cr contours. Both panels show the elongation of 
the contours, which imply some degeneracies between these parameters. 
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Fig. 3. — Probability spaces of a Leo show the degeneracy between stellar parameters. 
All the notations are the same as in the probability spaces of /3 Cas (see Fig. 2). The 
strong elongation of the contours in both panels suggest strong correlation between these 
parameters. The dashed lines show the vsmi range from McAlister et al. (2005), which 
intersects the probability contours with smaller common areas. 



-27- 



(3 Cas Image Reconstruction (3 Cas Model 




1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 

East (milliarcseconds) East (milliarcseconds) 



Fig. 4. — Images of /3 Cas . The left one shows the surface intensity distribution of 
/3 Cas from MACIM, overplotted with latitudes and longitudes from the model. The angular 
resolution is 0.57 mas(milli-arcsecond). The dashed contours represent the surface brightness 
temperatures of the image. The right one shows the image from model fitting, overplotted 
with brightness temperature contour from the model. The reduced of the images from 
MACIM and model fitting are 1.20 and 1.36. 
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Fig. 5. — Images of a Leo . The notations are all the same as those in images of /3 Cas (see 
Fig. 4). The angular resolution is 0.55 mas. The reduced of the images from MACIM 
and model fitting are 0.78 and 1.32. 
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Fig. 6. — P Cas positions on L-Rpr,] (left) and Hertzspr ung-Russell (right) diagrams based 
on model (lYi et al.ll200lU2003uDemarque et al.ll2004l ). In the left panel, the cross symbol 
with error bar stands for the rapidly rotating /3 Cas based on its Lboi and polar radius from 
modified von Zeipel model fitting. The square symbol with error bar is the non-rotating 
equi valent of Cas , t he corrections of Lboi and polar radius because of rotation is adopted 
from ISackmanru (119701 ). In the right panel, the solid line is the inclination curve, which shows 
how Lapp and T^pp change as a function of inclination angles. The star symbol is Cas with 
its estimated inclination angle. The meaning of the square symbol is the same as in the left 
panel (see Section 5). 



- 30 - 




Fig. 7. — a Leo position on L-Rpo\ and HR diagrams based on Y"^ model. The notations 
are the same as those on diagrams of /3 Cas (see Fig. 6). 
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Fig. 8. — The evolution of ste llar rotation. The model is adopted from the evolution 
of a 1.90 Mq non-rotating star (jPaxtonI |2004| . the web-browser interface is developed and 
maintained by Rich Townsend). The left panel: solid line is the ratio of the stellar radius 
to its value at the beginning of main sequence; dashed line is the ratio oi uj / ujait ((^ is 
angular velocity; Wcritis the critical angular velocity when the centrifugal force balances the 
gravity at the equator) to its initial value when the core and envelope are not coupled; dotted 
line is the ratio when they are totally coupled. The right panel: using the current u / ujcnt 
value (represented by asterisk) from model fitting, track back to its previous values assuming 
uncoupling and total coupling of the core and envelope. 
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Fig. 9. — Gravity darkening coefficient (/?) vs. temperature for four targets our group has 
studied. The sohd hne represents the theoretical relation between the gravity darkening 
coefficient and effecti ve temperature, adopted from the evolution of a non-rotating 2 solar 
mass star (jClaretll2000l ). The curve is extended to higher temperature for comparison with 
a Leo (see section 6.2 for details). The temperature range of each star contains temperature 
from the poles to equator. The range indicates the uncertainty from the model fitting. 
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Fig. 10. — Upper left panel: all seven nights visibility data of /3 Cas . The solid line and 
dotted line show the visibility curves of uniform disks with diameters of major and minor 
axis of /3 Cas from model fitting. The rest panels: the modified von Zeipel model (solid 
line) and MACIM image (dotted line, see section 4) vs. observed data (filled points with 
error bars) of /3 Cas from one single night. The reduced of model is 1.36 and that of 
image is 1.20. The eight data points in each sub-panel are from eight sub-channels of MIRC 
observation across H band. The x axis shows the wavelengths corresponding to the data 
points. The y axis shows which telescopes of CHARA have been used. 
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Fig. 11. — The similar panels of a Leo as those of /3 Cas in Fig. 10. The reduced 
model is 1.32 and that of image is 0.78. 



